REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 
Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection 
of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 

Final  Report 


1.  REPORT  DATE  (DD-MM-YYYY) 

02-06-2016 


4.  TITLE  AND  SUBTITLE 

Final  Report:  Modular  Universal  Scalable  Ion-trap  Quantum 
Computer 


3.  DATES  COVERED  (From  -  To) 

1 -Aug-20 10  -  31 -Jan-20 16 


5a.  CONTRACT  NUMBER 

W91  INF- 10- 1-0231 


5b.  GRANT  NUMBER 


6.  AUTHORS 
Jungsang  Kim 


5c.  PROGRAM  ELEMENT  NUMBER 
411359 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

Duke  University 

C/O  Office  of  Research  Support 

2200  W.  Main  St.,  Ste.  710 

Durham,  NC  27705  -4677 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS 
(ES) 

U.S.  Army  Research  Office 
P.O.Box  12211 

Research  Triangle  Park,  NC  27709-2211 


12.  DISTRIBUTION  AVA1LIBILITY  STATEMENT 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

58103-PH-MOC.l  12 


Approved  for  Public  Release;  Distribution  Unlimited 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

The  main  goal  of  the  original  MUSIQC  proposal  was  to  construct  and  demonstrate  a  modular  and  universally- 
expandable  ion  trap  quantum  computer.  This  architecture  has  two  separate  layers  of  scalability:  the  first  is  to 
increase  the  number  of  ion  qubits  in  a  single  trap,  with  full  control  over  the  qubits,  and  the  second  is  to  interconnect 
qubits  in  different  traps  using  photonic  channels.  The  proof-of-principle  demonstrations  for  single  and  two-qubit 
gates  in  a  single  ion  chain  had  been  carried  out  (on  a  2-ion  chain),  and  a  creation  of  entanglement  generation 

4fi  rt  4-  /%  ■**%  i  y-%  r»  *-%  *-%  1  r~i  iTrnr i  nl  n  y-~v  /I  *-%  n  4-  y-%  /I  \  \  7  y-%  1%  /-J  r»  «-%  •  -%-»  fi  lit  y-%  *  i  n  y>  y-~%  y-»  1  y-v  T  ■  *-%  4-  y-%  y>««n  4 «  *-%  y>  y-%  1  1  y-%  42  4-1-%  y-%  1-1  y-%  y-%  y-%  %-%— %  «-%  y-v  *-%  y-%  *-%  4n 


15.  SUBJECT  TERMS 

Ion  trap  quantum  computation,  scalable  modular  architectures 


—  i 


„U 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

15.  NUMBER 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF  PAGES 

UU 

UU 

UU 

UU 

Jungsang  Kim 


19b.  TELEPHONE  NUMBER 
919-684-3030 


Standard  Fonn  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39.18 


Report  Title 

Final  Report:  Modular  Universal  Scalable  Ion-trap  Quantum  Computer 

ABSTRACT 

The  main  goal  of  the  original  MUSIQC  proposal  was  to  construct  and  demonstrate  a  modular  and  universally-expandable  ion  trap  quantum 
computer.  This  architecture  has  two  separate  layers  of  scalability:  the  first  is  to  increase  the  number  of  ion  qubits  in  a  single  trap,  with  full 
control  over  the  qubits,  and  the  second  is  to  interconnect  qubits  in  different  traps  using  photonic  channels.  The  proof-of-principle 
demonstrations  for  single  and  two-qubit  gates  in  a  single  ion  chain  had  been  carried  out  (on  a  2-ion  chain),  and  a  creation  of  entanglement 
generation  through  photonic  channels  was  also  demonstrated.  We  had  an  ambitious  goal  of  integrating  all  of  these  components  and  assemble 
an  80-qubit  quantum  processor,  consisting  of  four  ion  chains  with  20  ions  in  each  chain,  interconnected  by  a  small  photonic  network.  While 
it  is  true  that  the  MUSIQC  collaboration  has  fallen  short  of  constructing  and  demonstrating  our  ambitious  goal  of  80-qubit,  networked 
quantum  processor  operation,  we  have  accomplished  significant  progress  in  the  field  of  ion  trap  quantum  computing  that  has  changed  the 
landscape  of  the  field  prior  to  MUSIQC  project.  We  have  developed  a  deep  understanding  of  fundamental  system  integration  challenges  that 
led  to  several  practical  solutions  we  have  been  able  to  partially  pursue  within  the  later  phases  of  the  MUSIQC  project,  as  well  as  new 
research  opportunities  outside  the  scope  of  the  MQCO  program. 
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Capstone  Report  for  IARPA  MQCO  Program 

Modular  Universal  Scalable  Ion-trap  Quantum  Computer  (MUSIQC) 

Summarized  by  Jungsang  Kim 
March  14th,  2016 

1.  Overall  Goals  for  the  MQCO  Capstone  Project 

The  goal  of  the  MQCO  Capstone  project  were  to 

(1)  Fully  characterize  the  error  contribution  for  our  prototype  systems  and  analyze  their 
impact  on  the  performance  of  the  algorithms 

(2)  Install  some  upgrades  to  our  existing  systems  to  improve  the  stability  of  the  system 
where  it  makes  sense 

(3)  Explore  a  range  of  advanced  quantum  control  techniques  to  improve  the  fidelities  of 
qubit  manipulation 

(4)  Improve  system  control  and  calibration  procedures  to  enhance  the  overall  error 
performance  of  our  algorithms 

(5)  Theoretical  activities  that  enable  design  and  implementation  of  adequate  quantum 
control  techniques  and  analyze  the  impact  of  the  fidelity  on  the  algorithm  execution. 

2.  Main  Accomplishments  in  the  Capstone  Project 

a.  University  of  Maryland  System 

At  University  of  Maryland,  we  established  the  experimental  capability  to  individually 
address  a  five-ion  chain  using  the  multi-channel  acousto-optic  modulator  (MC-AOM). 
The  gates  are  accomplished  by  breaking  the  control  Raman  pulses  into  nine  segments, 
which  is  designed  to  close  phase  space  loops  for  all  the  motional  modes.  The  fidelities  of 
the  Molmer-Sorensen  gates  among  arbitrary  pair  of  qubits  in  the  chain  were  characterized 
to  be  in  the  96-99%  range  (where  ~2%  SPAM  errors  were  subtracted),  as  shown  in 
Figure  1.  The  left  panel  in  Figure  2  shows  the  Rabi  oscillations  induced  on  each  ion  after 
the  chain  is  Doppler-cooled.  Due  to  the  counter-propagating  Raman  beam  geometry,  the 
Rabi  oscillations  are  sensitive  to  the  residual  ion  motion,  limiting  the  fidelity  of  the  single 
qubit  gates.  The  right  panel  in  Figure  2  shows  the  crosstalk  of  the  Raman  beams  to  the 
neighboring  ions.  The  crosstalk  to  the  nearest  neighbors  is  in  the  2-5%  range,  while  that 
to  the  next  nearest  neighbor  ions  is  below  1%  range.  There  are  three  sources  of  crosstalk: 
(1)  due  to  the  quality  of  the  optical  beam  at  the  ion  location,  (2)  due  to  RF  signal  leaking 
between  the  AOM  channels  and  (3)  due  to  the  acoustic  wave  leaking  between  the  AOM 
channels.  The  crosstalk  performance  can  be  improved  by  improving  the  design  and 
manufacturing  of  the  MC-AOM  device,  and  by  improving  the  shape  of  the  optical  beam 
incident  on  the  devices. 

Upon  completion  of  the  MS  gates,  we  proceeded  to  demonstrate  a  CNOT  gate 
between  the  all  possible  pairs  of  ions  in  the  chain.  This  is  realized  by  combining  the  MS 
gate  with  three  single-qubit  gates  applied  to  either  ion,  as  shown  in  Figure  3.  The  average 
time  for  the  gate  is  about  283  qs,  and  the  achieved  gate  fidelity  is  95.6%.  The  average 
detection  fidelity  is  about  98%  (SPAM  error  of  about  2%)  in  the  system. 


Figure  1:  M0lmer-S0rensen  gate  on  all  possible  ion  pairs  in  a  5-ion  chain  (less  the  left-right 
symmetry),  (a-f)  denote  the  two  qubit  gates  between  ion  pairs  1-2  (similar  to  4-5),  3-4  (2-3),  2-4, 
3-5  (1-3),  2-5  (1-4)  and  1-5,  respectively.  For  each  subfigure,  the  graph  on  the  left  shows  the 
population,  while  the  graph  on  the  right  shows  the  parity  oscillations.  The  gate  fidelity  is 
measured  to  be  in  96-99%  range,  when  the  SPAM  error  (~2%)  is  subtracted  out. 


Figure  2:  (Left)  The  Rabi  oscillations  of  each  qubit  in  the  chain,  when  the  chain  is  Doppler- 
cooled.  The  counter-propagating  gate  beams  make  the  Rabi  oscillations  sensitive  to  the  residual 
motion  of  the  ion,  limiting  the  fidelity  of  the  single  qubit  gates.  (Right)  Characterization  of  the 
crosstalk  in  Raman  beams.  The  crosstalk  between  the  nearest  neighbor  qubits  is  in  the  2-5%  range. 
While  that  for  the  next  nearest  neighbor  is  below  1%. 

The  transverse  mode  frequencies  of  trapped  ions  are  determined  largely  by  the  rf 
potential  delivered  to  the  trap  electrodes  through  an  rf  amplifier  and  a  helical  resonator. 
Both  amplifier  and  resonator  have  mechanical  and  thermal  fluctuations  that  affect  the  ion 
trap  frequencies,  and  are  a  leading  contributions  to  gate  infidelity.  We  have  stabilized  the 
secular  frequency  of  trapped  ion  motion  to  less  than  20  Hz,  at  a  trap  frequency  of  ~1 
MHz.  This  is  accomplished  by  capacitively  picking-off  1%  of  the  rf  high  voltage, 
rectifying  the  signal,  comparing  to  a  stable  voltage  source  acting  as  the  set  point,  with 
feedback  to  the  rf  source.  Both  pick-off  capacitor  and  diode  rectifier  circuits  (shown  in 
Figure  4)  are  designed  with  very  low  drift  and  temperature  sensitivity. 


(a) 


CNOT  [1:2]  F=96.4(6)%  CNOT  [1:3]  F=97.6(7)%  CNOT  [1:4]  F=95.9(7)%  CNOT  [1:5]  F=97.9(5)% 


(b) 


CNOT  [2:3]  F=95.6(6)% 


CNOT  [2:4]  F=98.4(7)% 


CNOT  [3:5]  F=97.6(6)% 


CNOT  [2:5]  F=96.8(7)% 


Figure  3:  Demonstration  of  a  CNOT  gate  between  the  all  possible  pairs  of  a  5-ion  chain,  (a) 
Schematic  of  the  CNOT  gate  realized  using  a  MS  gate  and  three  single-qubit  gates,  (b) 
Histograms  representing  the  population  of  the  input  and  output  states  for  the  CNOT  gate.  The 
gate  fidelities  range  from  95-99%  range,  without  the  contribution  from  SPAM  errors  (~2%). 


Figure  4:  (a)  Circuit  diagram  of  the  rectifier  that  processes  rf  pick-off  signal  from  the  helical 
resonator  for  feedback  and  comparison  to  set  point  potential,  including  dual  matched  diodes  for 
temperature-coefficient  compensation,  (b)  Mounting  of  rectifier  in  a  closed  SMA  package. 


Using  this  system,  we  have  implemented  the  Deutsch-Jozsa  algorithm,  the  Bemstein- 
Vazirani  algorithm,  and  the  quantum  Fourier  transform  algorithm  on  five  physical  qubits. 
Figure  5  shows  the  results  for  running  both  the  Deutsch-Jozsa  algorithm  (for  Figure  5  a  and 
b)  and  the  Bernstein- Vazirani  algorithm  (for  Figure  5c  and  d).  The  average  discrimination 
efficiency  is  0.967(2)  for  constant  and  0.932(3)  for  a  balanced  functions  in  the  Deutsch- 
Jozsa  algorithm.  For  the  four-qubit  Bernstein- Vazirani  algorithm,  the  single-shot 
detection  for  the  correct  estimated  sequence  is  achieved  with  an  experimental  average 
fidelity  of  0.903  (2). 


Balanced  Function 


Oracle 


c)  Ions  State 


Figure  5:  (a)  The  Deutsch-Jozsa  algorithm  on  5  ions.  Shaded  regions  show  reconfigurable  part  of 
the  sequence  where  the  light  and  dark  grey  sections  are  executed  in  separate  cases  to  achieve 
balanced  and  constant  functions  respectively.  CNOT  gate  combinations  (light  grey)  are  used  to 
create  parity  functions  of  upto  three  qubits.  A  rotation  on  X4  (dark  grey)  is  used  to  implement 
constant  function  f=l.  (b)  Measured  population  of  the  output  state  for  various  functions  post 
selected  on  X4  =  1.  Binary  outcome  of  the  state  111  indicates  that  the  function  is  constant,  while 
any  other  outcome  indicates  that  the  function  is  balanced,  (c)  The  Bemstein-Vazirani  algorithm 
circuit  where  shaded  region  contains  reconfigurable  CNOT  gate  combinations  used  to  implement 
different  oracle  states,  (d)  Measured  output  population  for  various  oracle  states.  The  output  is  the 
inverted  oracle. 


We  also  performed  a  fully  coherent  quantum  Fourier  transform  (QFT)  on  five  qubits 
using  techniques  that  can  easily  be  scaled  to  larger  registers.  We  examine  the 
performance  of  QFT  as  a  part  of  a  phase  estimation  protocol  without  implementing  the 
Oracle.  A  state  preparation  is  applied  instead  in  order  to  mimic  a  controlled  modular 
exponentiation  (Figure  6a).  The  state  preparation  involves  single  qubit  rotations  on  all  five 
qubits  which  is  followed  by  a  QFT  protocol.  Our  version  of  QFT  implementation  is 
different  than  earlier  trapped  ion  experiments,  in  the  sense  that  a)  we  perform  a  fully 
coherent  QFT  that  allows  us  to  implement  it  in  a  truly  reversible  way  as  opposed  to  a 
semi-classical  scheme,  b)  we  employ  two  qubit  gates  for  all  5C2  =  10  ion  pair 
combinations  in  the  register  thereby  illustrating  a  fully  connected  system  capable  of 
performing  any  arbitrary  gate,  c)  by  applying  composite  conditional  phase  gates  with 
varying  phases  for  all  pairs  of  ions,  we  demonstrate  the  modularity  and  hence  the 
scalability  of  this  scheme  for  implementing  the  QFT,  and  d)  by  avoiding  intennediate 
individual  detection  of  qubits  or  qubit  recycling,  as  required  in  a  semi-classical  scheme, 
our  coherent  implementation  is  free  of  accumulated  errors  due  to  imperfect  measurement 
and  decoherence  induced  by  the  measurement  process. 
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Figure  6:  Quantum  Fourier  transform  (QFT)  protocol,  (a)  Schematic  of  a  quantum  algorithm 
using  QFT.  Shaded  gray  region  represents  an  oracle  that  starts  with  a  superposition  state  applying 
Hadamard  (H)  gates  on  hypothetical  measurement  qubits  (dashed  red  wire)  followed  by 
controlled  unitary  (U)  operations  on  a  computational  quantum  bus  which  is  then  followed  by  QFT. 
Shaded  green  region  contains  experimental  sequence  where  solid  black  wires  represent  trapped 
ion  qubits.  ‘State  preparation’  contains  single  qubit  rotations  to  mimic  an  initial  state  given  by  an 
oracle  corresponding  to  period  finding/quantum  phase  estimation  algorithm.  ‘QFT’  contains  a 
sequence  of  conditional  phase  gates  and  rotations  (see  supplemental  material),  (b)  Quantum 
phase  estimation  using  five  measurement  qubits  is  performed  by  incrementally  setting  the  input 
phase  in  steps  of  2jt/64.  The  measured  output  is  a  state  whose  value  corresponds  to  the  input 
phase  is  detected  with  a  probability  >  0.6.  (c)  Quantum  period  finding  using  five  qubits.  Input 
states  are  prepared  using  single  qubit  rotations  to  create  amplitude  modulation  for  period 
{1,4,8,16,32}  and  phase  modulation  for  period  3.  The  squared  statistical  overlap)  SSO)  signifies 
the  fidelity  of  the  protocol  where  the  error  is  statistical. 


One  of  the  performance  limitations  on  the  application  of  multiple  gates  have  been 
identified.  The  4-photon  Stark  shift,  stemming  from  imperfect  polarization  of  the  355nm 
Raman  beams  and  driving  virtual  transitions  through  the  (1,1)  or  (1,-1)  Zeeman  levels,  is 
different  on  all  pairs  of  ions,  meaning  that  gates  on  particular  pairs  of  ions  must  account 
for  this  factor,  especially  when  qubits  are  involved  in  subsequent  gates.  We  have  finely 
optimized  the  Raman  beam  polarization  to  minimize  this  effect,  but  it  still  represents  a 
phase  of  ~2-5  degrees  on  each  gate,  presumably  from  vacuum  window  birefringence.  In 
the  future  this  is  not  a  problem  as  long  as  the  Stark  shift  does  not  drift  (and  it  appears  to 


be  stable  over  many  hours).  However,  we  currently  cannot  vary  the  phase  of  the 
individual  addressing  beams,  so  this  limits  the  number  of  gate  sequences  we  can  current 
apply.  We  have  implemented  a  5-qubit  fully  quantum  FFT  algorithm  on  a  variety  of 
input  states,  using  all  10  possible  2  qubit  gates  on  the  5  qubits.  We  are  currently 
analyzing  the  results. 


On  the  photonic  link  side,  the  Maryland  team  continues  to  press  on  the  establishment 
of  the  Ba-Yb  dual  species  system.  We  have  observed  (post-selected)  entanglement 
between  a  single  138Ba+  ion  and  the  polarization  of  an  emitted  photon  at  493nm  (see 
Figure  7a).  The  138Ba+  ion  is  weakly  (10%)  excited  with  a  500ns  pulse  of  resonant  light 
at  493nm,  and  we  collect  the  resulting  photon  through  an  NA0.6  objective,  orthogonal  to 
the  quantization  magnetic  field.  As  shown  in  Figure  7b  and  Figure  7c,  the  observed 
correlation  between  photonic  and  atomic  qubits  in  orthogonal  bases  proves  entanglement. 
Both  photonic  and  atomic  state  are  post-selected  (-10%  excitation  probability,  ~2% 
detection  efficiency  of  photon,  and  ~2%  detection  efficiency  of  1  lSBa+  qubit).  The  post- 
selected  fidelity  is  observed  to  be  F=0.86. 

138Ba+ 

Photon  Polarization  Rotation  138Ba+  Qubit  Rotation 


(a) 


(b) 
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Figure  7:  (a)  'NBa  energy  levels,  with  atomic  qubit  stored  in  the  electron  spin  states  in  the 
ground  state,  and  excitation  scheme  to  produce  single  photons  whose  polarization  (H  or  V) 
frequency  is  correlated  with  the  atomic  qubit  state,  (b)  Observed  correlation  between  atomic 
qubit  and  photon  polarization  in  two  bases  of  photon  measurement,  using  a  waveplate.  (c) 
observed  correlation  between  atomic  qubit  and  photon  polarization  in  two  bases  of  photon 
measurement,  using  a  waveplate. 


b.  Duke  University  System 

At  Duke  University,  we  pushed  on  two  fronts:  the  establishment  of  a  two-qubit  Molmer- 
Sorensen  (MS)  gate  on  one  of  our  chambers,  and  improved  photon  collection  into  a 
single  mode  fiber  using  the  high-NA  lens  for  better  ion-photon  entanglement  experiments. 

For  the  two-qubit  gate  experiment,  we  demonstrated  initial  two-qubit  MS  gate 
operation  using  a  single  pixel  PMT  detector,  with  histogram-fitting  to  distinguish  among 
the  cases  where  none,  single,  and  two  qubits  are  in  the  bright  state.  This  measurement 
method  works  for  the  two-ion  case,  but  not  in  the  multi-ion  case  we  are  aiming  for  in  our 
experiments.  Since  we  plan  to  get  to  longer  chains,  we  decided  to  install  a  multi-channel 
PMT  to  detect  the  qubit  state  of  each  ion  in  the  chain  independently.  For  this,  we  first 
realized  that  the  alignment  of  our  imaging  lens  to  the  ion  chain  was  very  far  from  optimal, 
leading  to  significant  aberrations  of  the  ion  image  that  spilled  photons  over  to  the 


neighboring  PMT  channels.  While  this  was  not  an  issue  in  a  single  PMT  experiment,  it 
accounted  for  about  9%  of  photon  spillage  to  the  neighboring  pixels,  and  about  3%  to  the 
next  nearest  neighboring  pixels  (which  we  used  to  detect  the  photons  from  the  adjacent 
ion).  We  went  through  a  realignment  of  the  imaging  optics  to  reduce  this  detection 
crosstalk.  Figure  8  shows  the  result  of  an  improved  alignment:  the  photon  spillage  to  next 
nearest  neighbor  was  reduced  to  below  1%.  A  careful  analysis  of  the  ion  image  still 
shows  a  small  halo,  leading  to  the  finite  spillage  to  neighboring  PMT  channels.  Since  we 
did  not  have  a  fine-tuning  capability  (of  motorized  stages,  etc.,  see  below),  we  could  not 
improve  the  alignment  better  than  this  limit.  The  photon  detection  crosstalk  among  the 
relevant  PMT  channels  will  lead  to  some  potential  crosstalk  error  in  qubit  detection  that 
will  depend  on  the  specifics  of  the  qubit  detection  protocol,  but  in  a  standard  histogram 
thresholding  approach,  we  expect  the  measurement  error  due  to  crosstalk  to  be  well 
below  1%  in  this  case. 


Channel 


Figure  8:  Measurement  of  pixel  crosstalk  for  the  multi-channel  PMT.  The  top  panel  shows  the 
photon  counts  from  each  channel  when  the  photons  collected  from  a  single  ion  is  aligned  to  the 
fourth  pixel  in  our  multi-channel  PMT.  While  there  is  significant  spillage  over  to  the  next  PMT 
channel,  the  spillage  to  the  next  nearest  channels  (pixel  2  and  pixel  6)  are  below  1%.  The  middle 
panel  shows  the  case  of  two  ions,  where  the  fluorescence  is  collected  to  channels  3  and  5, 
respectively.  The  inset  shows  the  image  of  ions  measured  by  the  EMCCD  camera,  where  a  small 
halo  around  each  ion  is  visible,  causing  the  spillage  to  neighboring  PMT  pixels.  The  bottom  panel 
shows  the  expected  photon  count  distribution  on  each  pixel,  based  on  the  measurement  of  the 
single  ion  case. 

Using  this  setup,  we  have  achieved  a  two  qubit  Molmer-Sorensen  gate  with  fidelities 
in  the  97%  range,  as  shown  in  Figure  9.  We  have  a  relatively  long  and  convoluted  Raman 
beam  path  from  the  355nm  Paladin  laser  to  the  ions,  where  one  of  the  beams  is  elevated 
about  50cm  above  the  optical  table,  and  above  our  DAC  electronics  (generating  quite  a 
bit  of  heat).  We  feel  it  is  difficult  to  keep  the  beams  stable,  and  therefore  the  Raman 
beam  intensities  at  the  ion  location  is  probably  fluctuating  at  a  few  %  level.  Our  MS  gate 


does  include  a  first-order  Walsh  correction,  but  it  is  not  very  effective  in  correcting  for 
intensity  errors.  We  have  made  efforts  to  stabilize  the  trap  frequency  as  the  detuning  of 
the  Raman  beams  also  drift  due  to  this  reason.  We  plan  to  implement  the  Maryland 
solution  on  trap  frequency  stabilization  in  the  future. 


Figure  9:  Characterization  of  a  two-qubit  Molmer-Sorensen  gate  at  Duke  University. 

Based  on  the  ray-tracing  simulations,  optimal  fiber  coupling  of  the  atomic  emission  at 
370nm  using  the  PhotonGear  lens  requires  critical  control  of  the  alignment  in  three 
degrees  of  freedom:  the  lens  position  along  the  optical  axis  (defined  to  be  the  z-axis,  to 
within  ~lpm),  and  x-  and  y-tilt  with  respect  to  the  optical  z-axis  (to  within  0.01°).  The 
alignment  is  completely  insensitive  to  rotation  around  the  z-axis  (cylindrical  symmetry), 
and  robust  against  x-  and  y-shifts  (~1 50um  object  space  over  which  the  image  is 
diffraction  limited).  The  true  z-axis  is  determined  as  the  normal  of  the  optical  window  for 
a  single  ion  case,  and  for  an  ion  chain,  we  require  that  the  chain  be  parallel  to  the  optical 
window  plane  to  within  <  1°  for  a  20  ion  chain.  The  initial  attempts  at  University  of 
Maryland  in  controlling  the  aberrations  have  been  very  effective  in  dramatically 
improving  the  coupling  efficiency,  but  it  looks  like  there  still  is  about  a  factor  of  3-4  to  be 
gained.  This  would  correspond  to  improving  the  entanglement  generation  rate  by  an  order 
of  magnitude,  so  it  is  definitely  worth  looking  into  further  improving  the  fiber  coupling 
efficiency. 

The  effort  at  Duke  started  with  an  adequate  design  of  motorized  translation  stage  and 
kinematic  mounts  that  provides  sufficient  resolution,  repeatability  and  stability  to 
position  the  PhotonGear  lens  at  the  optimal  location.  The  initial  alignment  process 
closely  followed  the  procedures  developed  at  Maryland,  but  with  an  improved  set  of 
alignment  optomechanics  and  diagnostic  probes  (high  resolution  imaging  and  comparison 
of  the  aberrated  images  with  ray  tracing  simulation  results).  After  the  final  correction,  the 
point  spread  function  (measured  as  the  shape  of  the  image  of  the  ion,  which  is  almost  an 
idea  point  source)  is  tightly  focused,  to  within  30%  of  the  diffraction-limited  spot  size. 


With  this  setup,  we  have  accomplished  coupling  50%  of  the  opolarized  light  collected 
within  the  NA  of  the  PhotonGear  lens  into  the  fiber  (see  Figure  10).  This  is  a  factor  of  2.5 
improvement  over  the  UMD  setup,  and  approaching  the  maximum  coupling  of  80% 
based  on  the  simulation  results  (determined  by  the  overlap  between  the  fiber  mode  and 
the  diffraction-limited  collection  mode  of  the  PhotonGear  lens).  This,  coupled  with  an 
superconducting  nanowire  single  photon  detectors  (SNSPDs)  that  have  been  tested  in  the 
Duke  lab  (this  part  of  the  project  is  being  funded  by  LPS/ARO)  which  shows  50%  system 
detection  efficiency  (-68%  device  efficiency  and  -73%  fiber  transmission),  we  anticipate 
achieving  the  ion-photon  entanglement  generation  rates  of  about  80  events  per  second,  a 
factor  of  17  improvement  from  the  recent  U.  Maryland  results.  The  next  step  is  to  verify 
the  fast  generation  rate  of  ion-photon  entanglement  generation.  In  order  to  facilitate 
efficient  diagnostics,  we  plan  to  use  the  polarization  qubit  for  the  photon,  which  is 
naturally  entangled  with  the  Zeeman  qubit  in  the  ion.  We  need  to  duplicate  this  setup  for 
the  second  chamber  before  ion-ion  entanglement  can  be  verified  at  higher  rates. 
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Figure  10:  Optimized  coupling  efficiency  of  the  emitted  photons  at  370nm  into  a  single  mode 
fiber,  measured  by  atomic  emission  driven  into  saturation,  (a)  The  detection  efficiency  of  the 
photons  scattered  by  174Yb+  ions  by  a  photomultiplier  tube  (PMT),  through  the  PhotonGear  lens. 
The  absolute  efficiency  can  be  extracted  from  the  saturation  curve  of  the  atomic  emission,  (b) 
Similar  curve  when  a  piece  of  single  mode  fiber  is  inserted  between  the  image  of  the  atomic 
emission  and  the  PMT.  Only  a- polarized  photons  couple  into  the  single  mode  fiber  in  this 
collection  geometry,  as  jr-polarized  photons  are  completely  rejected  by  the  mode  symmetry. 

c.  Georgia  Tech  Research  Institute  System 

The  main  goal  of  GTRI’s  effort  in  the  Capstone  project  was  to  push  through  major 
hardware  upgrades  in  the  system,  including:  (1)  Installation  of  a  300  mm  segment  of 
photonic  crystal  fiber  into  one  355  mn  Raman  gate  beam,  (2)  Reduce  the  size  of  that 
Raman  beam  into  the  AOM  used  to  turn  the  beam  on  and  off,  which  will  decrease  the 
pulse  rise  time  to  -  50  ns  (3x  improvement),  (3)  Redesign  path  for  the  other  Raman  gate 
beam,  such  that  it  is  no  longer  reflecting  off  a  dichroic  mirror,  (4)  Intensity  stabilization 
of  both  Raman  beams  with  high  bandwidth  feedback,  (5)  Study  of  our  feedforward  rep¬ 
rate  lock:  by  enclosing  the  entire  assembly  we  were  able  to  remove  spurious  noise  peaks, 
and  (6)  Installation  of  low  noise  current  drivers  for  magnetic  field  coils. 


Within  the  first  two  months,  we  completed  the  installation  of  the  fiber  segment, 
reduced  the  Raman  beam  into  the  AOM,  fixed  the  beam  paths,  and  implemented  the 
power  stabilization  feedback  loop  for  both  Raman  beams.  Upon  our  initial  single  qubit 
Rabi  oscillations  experiment,  we  realized  that  the  high  power  used  in  the  beams  have 
introduced  additional  qubit  frequency  shift  due  to  the  self-resonant  Raman  transition 
16MHz  from  the  carrier.  We  modeled  this  shift  and  determined  that  the  PB2  pulse 
sequence  removes  most  of  the  error  ( 1  O'5)  on  targeted  ions.  However  untargeted  ions 
experience  a  measurable  Rz  rotation.  We’ve  incorporated  this  knowledge  into  our  gate 
compiler,  such  that  it  applies  the  appropriate  rotations  on  the  whole  ion  chain.  Our 
simulator  also  accounts  for  the  unwanted  Rz  rotations.  Having  accounted  for  these  errors, 
we  utilized  Randomized  Benchmarking  to  characterize  single-qubit  gate  fidelities.  The 
infidelity  of  a  single-qubit  gate  operation,  in  a  two-ion  chain,  is  still  around  two  percent. 
To  determine  what  may  be  causing  this  error,  we  utilized  the  “quantum  lock-in” 
technique  with  80  echo  n  pulses.  These  scans  indicate  a  large  noise  source  near  45  kHz. 

We  spent  a  lot  of  time  tracking  down  the  source  of  this  noise.  In  the  process,  we 
discovered  that  our  repetition  rate  lock  for  the  355nm  laser  had  some  residual  phase  noise 
in  ~lms  timescales,  which  was  also  observed  at  Duke  and  Sandia  in  the  past.  We  tracked 
down  that  and  fixed  the  problem,  so  that  the  coherence  time  measurement  (Ramsey 
spectroscopy)  was  recovered  to  expected  levels.  However,  this  did  not  resolve  the  45kHz 
noise  issue.  We  have  a  pretty  good  understanding  of  what  this  might  be:  a  two-phonon 
transition  in  the  ion  chain,  where  the  ion  absorbs  (emits)  into  the  rocking  mode  and  emits 
(absorbs)  into  the  COM  mode.  As  the  two  radial  modes  are  separated  by  only  45  kHz,  the 
net  transition  is  45  kHz  detuned  from  the  carrier.  While  it  is  a  second  order  coupling,  our 
62.5  kHz  Rabi  rate  is  larger  than  the  detuning.  This  error  can  be  decoupled  by  modifying 
the  pulse  sequences,  or  mitigated  with  better  cooling. 

Although  we  have  not  been  able  to  resolve  all  of  these  experimental  issues  that  were 
uncovered  during  the  Capstone  period,  we  have  carried  out  the  3-qubit  Bemstein- 
Vazirani  algorithm  at  the  end  of  the  Capstone  period.  The  average  success  probability  of 
the  algorithm  has  generally  improved  greatly  over  the  experiments  at  the  end  of  the  5- 
year  MQCO  program,  as  shown  in  Figure  1 1 . 


Figure  ll:The  success  probability  of  3-qubit  Bernstein-Vazirani  algorithm  in  GTRI  system.  The 
performance  for  each  oracle  state  showed  significant  improvements  over  previous  experiments 
performed  with  a  setup  before  the  Capstone  upgrades. 


d.  Theory  Efforts  in  Noise  Modeling 

The  challenge  of  disentangling  the  spin  and  motion  degree  of  freedom  is  crucial  in 
implementing  a  high  performance  gate.  We  developed  theoretical  model  allowing  spin- 
motional  disentanglement  to  be  quantified  in  the  presence  of  phase-modulated  MS 
gates.  Theory  provides  a  good  prediction  of  expected  MS  lineshape  as  a  function  of 
Raman  detuning  from  motional  modes.  We  have  experimentally  verified  that  the  phase- 
shifted  gate  design  works  well,  and  started  looking  at  its  application  to  ion-motion 
entanglement.  However,  we  are  seeing  that  there  is  residual  micromotion  that  cannot  be 
corrected  in  our  setup,  making  it  difficult  to  simultaneously  cool  both  the  red  and  the  blue 
sidebands  of  noise.  This  could  be  arising  from  axial  micromotion,  potentially  due  to  a 
slight  misalignment  of  the  two  end  caps  in  our  macroscopic  trap,  and  the  lack  of  in-situ 
RF  filtering  (the  end  cap  electrodes  picking  up  a  small  portion  of  the  RF  from  the  RF 
electrodes). 

We  have  recently  obtained  a  better  estimation  of  phase  noise  of  the  oscillators 
generating  the  gate  pulses  on  the  qubit  gate  perfonnance,  and  conclude  that  the  impact  of 
the  phase  noise  could  be  much  larger  than  expected,  when  controllers  with  wide 
bandwidths  are  desired.  This  is  because  the  phase-noise  far  from  the  carrier  has  a  much 
higher  contribution  to  the  qubit  noise  than  previously  considered:  making  the  controller 
narrow-band  helps  with  the  noise  situation,  at  the  expense  of  the  control  bandwidth. 

Randomized  benchmarking  widely  used  for  characterizing  the  error  probabilities  of 
quantum  gates  lets  you  estimate  the  average  error  rates,  while  the  fault-tolerance 
threshold  theorems  require  a  stricter  notion  of  worst-case  error,  measured  by  the  diamond 
distance.  The  gap  between  the  average  error  and  worst-case  error  could  be  very  large:  in 
the  worst  case  of  unitary  (or  systematic)  errors,  this  gap  could  be  many  orders  of 
magnitude  in  the  low  error  regime.  We  have  developed  a  theoretical  framework  that 
allows  one  to  clearly  distinguish  the  extent  of  unitary  errors  present  in  the  system,  with  a 
goal  of  identifying  and  eliminating  them.  In  this  limit  where  all  coherent  errors  are 
eliminated,  the  average  and  worst-case  errors  now  coincide,  and  the  error  probabilities 
estimated  by  randomized  benchmarking  is  directly  relevant  for  the  realization  of  fault- 
tolerant  error  correction.  At  a  more  concrete  level,  we  studied  the  effect  of  modeling 
errors  with  extensions  to  the  Pauli  error  model  and  how  the  choice  of  approximation 
changes  the  behavior.  For  incoherent  errors  we  have  found  Pauli  channels  optimized  on 
fidelity  are  good  models  even  for  non-unital  errors  such  as  amplitude  damping  in  the 
context  of  error  correcting  codes.  A  purely  coherent  error  model  of  unwanted  rotation  by 
a  fixed  angle  after  each  gate  is  not  well  captured  by  either  Pauli  model.  The  Pauli 
channel  optimized  on  fidelity  (a  Pauli  Twirling  Model)  results  in  an  optimistic  threshold 
while  the  Pauli  channel  optimized  with  the  constraint  that  the  distance  for  pure  states 
must  be  worse  than  for  the  real  channel  results  in  a  pessimistic  threshold.  This  can  be 
understood  by  looking  at  the  leading  order  term  in  the  distance  for  small  distances. 

e.  Theory  Efforts  in  Algorithms  Development  and  Novel  Architectures 

We  (El.  Michigan)  carried  out  more  detailed  simulation  for  realization  of  the  boson 
sampling  algorithm  with  trapped  ions.  A  moderate  system  size  with  20  to  30  ions 
provides  the  most  interesting  physical  region  for  this  test  to  beat  the  classical  computers 
on  this  specific  problem. 


The  coherent  and  dissipative  processes  (such  as  photon  entanglement  or  qubit  state 
initialization  and  detection)  that  must  coexist  in  a  modular  architecture  like  MUSIQC 
calls  for  a  second  ion  species,  either  to  support  sympathetic  cooling  or  to  establish  the 
remote  entanglement  generation  using  photons  at  a  different  wavelength.  We  propose  a 
new  way  for  spectral  addressing  in  an  ion  crystal.  In  order  to  avoid  resonant  photons  to 
scatter  and  decohere  the  qubit  ions,  one  approach  is  to  have  the  ancilla  ions  represented 
by  a  different  ion  species,  and  such  as  a  dual-ion  species  system  the  MUSIQC  team  has 
been  focusing  on.  Apart  from  the  laser  complexity  to  control  two  species  of  ions,  an  issue 
for  this  approach  is  that  it  is  difficult  to  generate  a  regular  distribution  pattern  of  the 
ancilla  ions  around  the  computational  ions,  and  such  a  regular  pattern  is  likely  required 
for  realization  of  fault- tolerant  quantum  computation.  We  propose  a  different  approach  to 
address  this  problem,  where  the  ancilla  ions  are  represented  by  the  same  species  of  ions, 
but  on  a  different  metastable  state  (the  shelter  D  level).  We  apply  laser  cooling,  repetitive 
readout  measurements,  and  probabilistic  entangling  operations  all  on  these  D  level  ions, 
which  leave  the  laser  beams  and  spontaneous  emission  photons  having  no  influence  on 
the  computational  ions.  It  is  easy  to  get  a  regular  pattern  of  these  ancilla  ions  through 
spatially  focused  laser  beams,  which  bring  the  ancilla  ions  to  the  D  level. 

We  also  studied  the  possibility  of  constructing  robust  distributed  entanglement 
between  small  ELUs  similar  to  those  considered  in  the  MUSIQC  program,  that  might  be 
suited  for  a  long-term  application  to  communications,  with  the  following  properties: 

1.  With  L  the  total  distance  between  the  communicating  parties,  the  operational  cost 
is  of  order  L  log2  (L). 

2.  Under  the  assumption  of  fast  classical  processing  (one-way  communication  + 
local  computation),  the  rate  of  secure  communication  is  independent  of  L. 

3.  The  communication  structure  can  be  conveniently  built  from  5 -qubit  ELUs, 
similar  to  those  considered  in  MUSIQC  scalability  challenge. 

Compared  to  the  recent  proposals  on  QECC-based  third-generation  quantum 
repeaters,  the  present  construction  provides  the  compatibility  with  small  ELUs. 
Traditionally,  it  is  believed  that  cramming  all  the  necessary  components  into  small 
repeater  stations  has  a  tendency  to  result  in  loss  of  performance:  we  find  that  this 
degradation  is  not  necessary.  Similar  to  the  solution  of  MUSIQC  scalability  challenge, 
we  can  utilize  3D  cluster  states  to  get  the  above  3  characteristics.  In  this  construction,  we 
(1)  build  a  3D  cluster  state  stretching  from  one  communication  party  to  the  other,  using 
5-qubit  ELUs,  as  shown  in  Figure  12.  Four  of  the  qubits  in  each  ELU  are  used  to 
establish  Bell  pairs  with  neighboring  ELUs  via  photonic  links.  (2)  Locally  measure  all 
qubits  in  the  3D  cluster  in  the  eigenbasis  of  ox,  except  for  one  at  the  respective  locations 
of  each  communicating  party.  (3)  From  all  the  measurement  outcomes  construct  the 
error-syndrome  and  correct  for  the  most  likely  compatible  error.  This  produces  an  almost 
pure  Bell  state  between  the  parties. 

The  protocol  works,  of  course,  only  below  the  fault-tolerance  threshold  for  3D  cluster 
states,  which  has  numerically  been  established,  for  a  large  range  of  error  models,  to  be 
(Ka)  »  0.72,  where  Ka  is  an  elementary  cluster  state  stabilizer  operator.  From  this 
criterion,  we  obtain  the  following  bound  on  the  error  rates, 


Therein,  s  is  the  gate  error  (all  one-and  two-qubit  gates),  T  is  duration  of  slowest  gate  and 
rD  is  the  decoherence  time. 


3D  cluster  state 


Figure  12:  3D  cluster  state  for  fault-tolerant  entanglement  distribution  and  quantum 
communication. 


